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Abstract. We show that the quantum coherent transfer of excitations between 
biomolecular chromophores is strongly influenced by spatial correlations of the 
environmental fluctuations. The latter are due either to propagating environmental 
modes or to local fluctuations with a finite localization length. A simple toy model 
of a single donor-acceptor pair with spatially separated chromophore sites allows 
to investigate the influence of these spatial correlations on the quantum coherent 
excitation transfer. The sound velocity of the solvent determines the wave lengths 
of the environmental modes, which, in turn, has to be compared to the spatial 
distance of the chromophore sites. When the wave length exceeds the distance 
between donor and acceptor site, we find strong suppression of decoherence. In 
addition, we consider two spatially separated donor- acceptor pairs under the influence 
of propagating environmental modes. Depending on their wave lengths fixed by the 
sound velocity of the solvent material, the spatial range of correlations may extend over 
typical interpair distances, which can lead to an increase of the decohering influence 
of the solvent. Surprisingly, this effect is counteracted by increasing temperature. 
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1. Introduction 



The photosynthetic conversion of physical energy of sunlight into its chemical form 
suitable for cellular processes involves many physical and chemical mechanisms [H E]. 
Photosynthesis starts with the absorption of a photon by a light-harvesting pigment 
forming an exciton, followed by the transfer of the exciton to the reaction center, where 
charge separation is initiated. It is the nature of this transfer process which is presently 
in the focus of intense research. Recent experiments [Sj H] provided evidence that an 
incoherent hopping model seems not to be sufficient to describe long-lasting beating 
signals in a two-dimensional Fourier transform electronic spectrum [5] recorded from 
green sulfur bacteria such as Chlorobium tepidum. Here, the energy transfer between the 
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main chlorosome antenna and the reaction centers is mediated by the Fenna-Matthews- 
Olson (FMO) protein O [71 E] which contains bacteriochlorophyll (BChl) molecules. 
The FMO protein is a trimer made of identical subunits, each of which contains seven 
BChl molecules and no carotenoids. Due to its small size, it represents an important 
model system for photosynthetic energy transfer and has been extensively studied 
experimentally and theoretically. The observed |3] long lived electronic coherence lasted 
up to time scales comparable to the time scale of the energy transport. The experiments 
were performed at low temperature T = 77 K and clearly suggest that the exciton moves 
coherently through the FMO complex rather than by incoherent hopping. Similarly, Lee 
et al. |1] found coherent beating signals at low temperatures in a two-color electronic 
coherence photon echo experiment. It allows to directly probe electronic coherences by 
mixing of the bacteriopheophytin and accessory bacteriochlorophyll excited states in the 
reaction center of the purple bacterium Rhodobacter sphaeroides. These measurements 
were performed at 77 K and at 180 K. The coherence beatings in these measurements 
can only be explained [9] by a strong correlation between protein-induced fluctuations 
in the transition energies of neighboring chromophores, leading to the conclusion that 
protein-correlated environments in fact preserve and support electronic coherence in 
photosynthetic complexes. 

Recently, an ultrafast polarization experiment [10] has revealed quantum coherent 
intrachain (but not interchain) electronic energy transfer in conjugated polymers with 
different chain conformations as model multichromophoric systems at room temperature. 
The data suggest that chemical donor-acceptor bonds help to correlate dephasing 
perturbations. By introducing an angle-resolved coherent optical wave mixing technique, 
the quantum beating signals between coherently coupled electronic transitions in the 
light-harvesting complex of purple bacteria are directly observed [llj. Also in the 
light-harvesting proteins of cryptophyte marine algae, quantum coherent couplings have 
been identified by exceptionally long-lasting excitation oscillations even at ambient 
temperature [T2] . 

The consequences of these seminal experiments on long-lived electronic coherences 
are immediate: the excitation can move rapidly and reversibly in space, allowing 
for a very efficient search for an energetic minimum in photosynthesis. The 
subsequent trapping of the excitation, however, must be optimized by properly 
adjusted environmental fluctuations. In view of potential applications, the recent 
progress helps to understand the design principle of photosynthetic complexes [13] 
and to exploit the near-unity efficiency of energy transfer which is believed to result 
from the constructive interplay of quantum coherence and slow, spatially correlated 
environmental fluctuations. This could open the door to efficient future artificial 
light-harvesting complexes flnding applications in optimized organic solar cells. The 
combination of optimized exciton trapping [H] with powerful quantum coherent 
adaptive control schemes [151 HSl E] could, in addition, allow to exploit quantum 
effects to direct the outcome of photochemical processes. Effects beyond the lowest- 
order Forster treatment, such as the failure of the point-dipole approximation and the 
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ensuing solvent screening and the sharing of common bath modes have been considered 
in a generahzed Forster theory [18]. 

Theoretically, photosynthetic energy transfer processes in light-harvesting 
complexes are often discussed in simplified low- dimensional models describing a few 
individual chromophores which mutually interact by dipolar couplings and which 
are exposed to the fluctuations of the polar solvent molecules and the protein host 
[H |2l [191 IMl HH [221 [23]. Two limiting cases are commonly considered: (i) When 
the dipole coupling between the chromophores is weak in comparison to the coupling 
of the chromophore to environmental fluctuations, the excitons are considered to 
be localized at the chromophore sites. The weak electronic coupling can then be 
treated perturbatively, resulting in an incoherent hopping dynamics described by the 
standard Forster theory [24]. (ii) In the opposite limit of weak coupling to the 
environmental fluctuations, standard perturbative quantum master equations are used 
[25l[26l[2ll[2Hl[29l[30l[ni[33[33l[33[3g resulting in a damped coherent dynamics 
for the exciton transfer. They naturally are based on the Markovian approximation, 
which renders the time evolution memoryless and allows for a straightforward numerical 
solution. These types of master equations are appropriate in the case of a clear 
separation of time scales, i.e., when the time scale on which the environmental 
fluctuations occur is much smaller than that on which the system dynamics evolves. 
Formally, this is captured by the requirement that the bath reorganization energy 
[20] is much larger than the typical system energy. This condition is typically not 
fulfilled for the energy transfer dynamics in biomolecular light-harvesting complexes 
in a protein-solvent environment [2T| [22| [23| [36] . Here, in addition to the rather 
slow polarization fluctuations of the polar solvent molecules, the protein cage acts 
also as a frequency filter which particularly shapes the frequency distribution of the 
environmental modes. A similar effect arises for electronic spin qubits in semiconductor 
quantum dots [37] and in donor-based charge qubit crystal systems [38]. Here, the 
designed geometrical shape of the system structures the acoustic phonon spectrum and 
gives rise to strong non-Markovian effects. In the context of excitonic energy transfer, it 
has recently been shown [391 [101 [S] that under realistic biomolecular circumstances, 
time-local master equations become increasingly unreliable when both time scales 
become comparable. In addition to the violation of the Markovian assumption, the 
coupling between chromophores and environment cannot be considered as weak enough 
to allow for a lowest order perturbative treatment [211 [221 [23]. It was recently shown 
that the latter condition renders any weak coupling approach questionable [H]. Non- 
Markovian approaches have been employed beyond a lowest-order treatment [39] and by 
coupling each chromophore to a single damped harmonic mode [35], whose existence, 
however, was not further motivated. Alternatively, numerically exact simulations of 
the real-time dynamics of quantum coherent energy transfer under realistic conditions 
have been carried out [IHl [HI [12] by employing the quasiadiabatic propagator path- 
integral (QUAPI) [131 HH HH]- By this, it has been shown that the rather slow 
polarization fluctuations are one possibility to enhance quantum coherence in the 
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transfer processes. The coupling of two chromophore pairs to the common slowly 
fluctuating modes even allows to entangle two excitonic pairs over surprisingly long 
times even at room temperature |lOj. However, when the fluctuations are fast, no 
entanglement is created even when the two pairs couple to the same modes. The 
possibility of entanglement and the role of non-Markovian contributions in biomolecular 
complexes have also been re-addressed in recent works |33l IMl US]. Slow fluctuations 
have also been treated |17] by different variants of cumulant expansion techniques 
and by statistical averaging approaches over static disorder in sum-over-eigenstates 
approaches. Recently, the standard Redfield equations, which are valid in the weak- 
coupling regime, have been extended by generalizing the Redfield relaxation tensor on 
the basis of the Lindblad quantum master equation [58]. This technique goes beyond 
the secular approximation and thus can include effects of stronger coupling. However, 
the approach is still memoryless and leads to time-local evolution equations. 

Despite the fact that the experimental coherence beatings [9] could only be 
explained including strong correlations between protein-induced fluctuations in the 
transition energies of neighboring chromophores, the influence of these correlations 
received little attention in the theoretical investigations. Nazir |19] investigated the 
influence of correlated fluctuations on a donor-acceptor system for strong system bath 
coupling. Correlations in a super-Ohmic bath are found to suppress the crossover 
to incoherent dynamics at high temperatures, which is in line with the common 
expectation [20] that a super-Ohmic bath naturally provides only weakened influence of 
the fluctuations on the system. Similar effects would be expected for correlations in the 
environment of the chromophores whose spectrum is typically assumed to be Ohmic. 
Fassiolo et al. [29] discussed the influence of correlations on the trapping probability in 
a ring of chromophores within a Lindblad master equation approach. They find that the 
correlations between environmental fluctuations allows to tune the trapping probability. 

From a condensed matter point of view, environmental fluctuations are normal 
modes of the bulk material and thus are typically phonons which propagate through 
the material. Accordingly they couple to all chromophores with amplitude differences 
determined by the phase differences due to finite times the modes need to propagate from 
one chromophore site to the next. However, rattling of side chains of macromolecules 
in these highly disordered protein environments might well be viewed as a localized 
excitation. In most existing studies it is assumed that each site of a multichromophoric 
array is coupled to its local environment. Including in these approaches spatial 
correlations could be achieved by assuming a finite localization length of the excitations. 

In the next section we discuss the influence of spatially correlated environmental 
fluctuations on a donor-acceptor model using the numerical exact quasi-adiabatic path 
integral propagator approach which allows us to treat realistic strong couplings and slow 
environments. Although being a clear oversimplification to realistic exciton transport 
the donor-acceptor model serves as a toy model to study the influence of spatial 
correlations and the difference between propagating and localized modes on quantum 
coherence in detail. In the third section we will discuss how spatial correlations influence 
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two donor-acceptor pairs which are initially uncoupled. Depending on the distance r^a 
between donor and acceptor and the distance r between the two donor-acceptor pairs 
spatially correlated fluctuations increase or decrease the decay rates of the coherent 
dynamics. Finally we discuss and summarize our results. 

2. Spatial environmental correlations in a single chromophore pair 
2.1. Model 

The simplest way to model a single chromophore (or pigment) is by describing it 
as a quantum two-level system consisting of a ground and an excited state, which 
are separated by the energy gap e. When the electron is in the excited state, it 
is localized by its attractive interaction with the hole it left. This dipole electron- 
hole configuration forms an exciton. A formal description can be given in terms of 
the Pauli matrix r^. Environmental fluctuations will cause transitions between the 
ground and the excited state and will add a fluctuating energy. Experimentally it 
is known that the recombination time is of the order of nanoseconds, whereas the 
complete energy transfer through the complex is of the order of picoseconds. Thus 
the environmental fluctuations causing recombination are negligible. Describing the 
fluctuations by harmonic oscillators, which couple linearly to the chromophore, results 
in the independent boson model for a single chromophore 

H = e| -|e)(e|5^A,(r)g, + W 

k k 

where we introduced the position and momentum operators, and pk, of the mode 
with wave vector k and its coupling Ak(r) to the chromophore which depends on the 
amplitude of the fluctuation at the position r of the chromophore. We explicitely coupled 
the environmental fluctuations only to the excited state, |e), which expresses the fact 
that the electronic ground state energy is defined by including all vibrational equilibrium 
energies. We fixed h = = ^ which we keep below. 

We are not aiming at a complete description of exciton transfer dynamics in 
complexes like FMO but are merely interested in the question how spatially correlated 
environmental fluctuations influence the transfer process between two excitonic sites. 
For the sake of simplicity, we restrict the model under consideration to two chromophore 
sites (acceptor and donor at Ya/d) which contain a single exciton. We do not consider 
different site energies (i.e., = e^) and are thus lead to the donor-acceptor Hamiltonian 

i/da = ^A{|rf)(a| +\a){d\} + Y, K)(^lE^k(r.)gk+^$^(p^ + c^M (2) 

i=a/d k k 

of a single chromophore pair. The state \d) (|a)) denotes the exciton to be at the donor 
(acceptor) and A is the respective dipole coupling matrix element |21] . 

Comparison of the donor- acceptor model with the spin-boson model 
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The donor-acceptor Hamiltonian, Eq. ([2]), can easily be transformed into a Hamiltonian 
which is closer to the widely studied spin-boson model 120] 

k k 

+ ^ E + ^k?k) (3) 

k 

by introducing the Pauh matrices {l,ax,ay,az} with = \d){a\ + \a){d\ and = 
\d){d\ — \a){a\. 

One difference is given by the term proportional to Ak(ri) + Ak(r2) which couples 
to the identity operator 11 of the donor-acceptor system thus causing fluctuations of 
the reference energy of the donor-acceptor system. Accordingly it is irrelevant for 
its dynamics. However, this term modifies the bath modes by shifting their zero- 
point energies and thus changes their thermal equilibrium state. In the spirit of 
dissipative quantum dynamics, the treatment of system-bath problems typically rely 
on the assumption that the bath is only weakly influenced by the coupling to the 
system itself and thus the mentioned effects should not affect the dissipative dynamics 
of the donor-acceptor system. However, at strong coupling or for a slow bath with 
cut-off frequency Uc ^ A the validity of these assumptions is questionable [50]. Thus 
the dynamics generated by our donor-acceptor Hamiltonian in Eq. ([2]) differs from the 
standard spin-boson problem in three aspects. First, the thermal equilibrium state 
to which the total system is driven is different and, second, the factorized initial 
conditions for both cases reflect two different initial conditions. Third, in the spin- 
boson model, a single bath is coupled to the system whereas in the donor-acceptor 
Hamiltonian, two baths are coupled to the system states. Thus, even when both baths 
are mutually uncorrelated, the resulting rates for the donor-acceptor system are twice 
as large compared to the spin-boson model (assuming that all baths are coupled with 
equal strength). This has to be taken into account when comparing results from both 
approaches. 



2.2. The quasiadiabatic propagator path integral for the multi-bath case 

The dynamics of the donor-acceptor is characterized by the time evolution of the reduced 
density matrix p{t), which is obtained after tracing out the environmental (or bath) 
degrees of freedom, i.e., 

p(t) = Tr{f/(t,0)iy(0)f/-i(t,0)}^ (4) 

and 

f/(t,0) = Texpj--^^ dsHda^ . (5) 

Here, U (t, 0) denotes the propagator of the full system plus bath and T denotes the 
time-ordering operator. W{0) is the total density operator at initial time set at t = 0. 
We assume standard factorizing initial conditions [20], i-e. 1^(0) oc p(0) exp(— if^/T), 
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where the bath with the Hamiltonian Hb = \ 'Ylik ^1 + '^llD thermal equihbrium 
at temperature T and the system is prepared according to p(0). Throughout this work, 
we always start with the exciton at the donor site, i.e., p(0) = \d){d\. 

We calculate p(t) using the numerically exact quasiadiabatic propagator path- 
integral (QUAPI) [ini HH SH] scheme. For details of the iterative technique, we refer 
to previous works [13], HU HS] . In brief, the algorithm is based on a symmetric Trotter 
splitting of the short-time propagator /C(tfc_(_i, tjt) for the full Hamiltonian into a part 
depending on the system Hamiltonian and a part involving the bath and the coupling 
term. The short-time propagator describes time evolution over a Trotter time slice 
6t. This splitting is by construction exact in the limit 5t — t- but introduces a finite 
Trotter error for a finite time increment, which has to be eliminated by choosing 6t 
small enough such that convergence is achieved. On the other side, the bath degrees 
of freedom generate correlations being non-local in time. For any finite temperature, 
these correlations decay exponentially fast at asymptotic times, thereby setting the 
associated memory time scale. QUAPI now defines an object called the reduced density 
tensor, which lives on this memory time window and establishes an iteration scheme in 
order to extract the time evolution of this object. Within the memory time window, 
all correlations are included exactly over the finite memory time r^em = K5t and can 
safely be neglected for times beyond Tmcm- Then, the memory parameter K has to be 
increased, until convergence is found. The two strategies to achieve convergence are 
naturally countercurrent, but nevertheless convergent results can be obtained in a wide 
range of parameters. 

For the purpose of this work, we have to extend the standard formulation of QUAPI 
which only includes the coupling to one bath. The entire influence of a single bath 
coupled via the operator s to the donor-acceptor system is described in terms of the 
real-time path-integral formulation by the influence functional 

I{{st, sT}- 5t) = exp j-i ^ [st - S-] - r,t,s^] \ (6) 

where the path segments sf associated to a Trotter time slice i given as interval 
[{i -\)5t,{i+ \)5t\ (with total time t = N6t) are assumed to have constant values 
over a single time slice. The number of path segments within a Trotter time slice is 
given by the dimension of the Hilbert space in which the system-bath coupling operator 
lives. The superscript ± denotes the propagation direction forward or backward in time 
since we work with density operators. The total path integration over all paths s^(t') has 
to be performed as the discrete sum over all configurations {sf,s^} of paths segments 
between initial and final time. The time-discrete bath correlators rjai are defined in 
Ref. [13] and the superscript * denotes the complex conjugate. 

Multiple independent baths, HBa, which couple to system operators Sa will simply 
cause a product of influence functionals since each bath acts separately as described 



Quantum coherent hiomolecular energy transfer with spatially correlated fluctuations 8 
above. Thus, the total influence functional assumes the form 



Here we denoted the bath correlators 7]l^, with the additional superscripts since for 
differing baths the correlators will differ. 

The question whether the environmental fluctuations act locally or in a correlated 
manner can be tackled by the following extension. Local fluctuations couple to the donor 
and to the acceptor separately and independently. This implies that Eq. ([7]) describes 
all effects due to environmental fluctuations. If, however, the fluctuations are caused by 
extended waves, like phonon modes or if the fluctuations rattling the donor can at least 
partially still be felt at the acceptor site, then the fluctuations at the various sites are 
no longer independent and spatial correlations have to be taken into account. Hence, 
Eq. ([7]) has to be generalized to 



where 77^^ are the mixed bath correlators expressing the correlations of the fluctuations 
acting at operator s^ and s^. The detailed numerical evaluation of the influence 
functional in the extended QUAPI scheme becomes more involved, but the general 
procedure is not affected by this extension. 

2. 3. Correlated environmental fluctuations at different sites 

In a crystal, environmental fluctuations acting, e.g., on electrons, are generated by 
vibrations of the lattice atoms and are the well-known phonons. Phonons are also 
present in disordered media (condensed, soft or fluid). In the sense of propagating 
modes of the host material which evolve with time through the medium, they are 
commonly limited to the low energy sector or, more specifically, to energies associated 
to wave lengths on which the material appears homogeneous. Once the wave length 
becomes smaller than the disorder length scale, the modes can generally be thought of 
as localized fluctuations. However, their localization length (or radius) is still connected 
to the wave length of the mode and accordingly even a localized mode extends over some 
finite spatial range. The same picture holds for dipolar fluctuations in solvents and 
vibrations of charged macromolecular side chains forming the bio-environment of light 
harvesting complexes. Excitons, as coupled electron hole pairs, have an electric dipole 
moment fl which couples to the electric field E{r) at the exciton position r generated 
by the environmental dipolar fluctuations (as illustrated in Fig. [1]). This results in the 
interaction energy Hsb = pl ■ E{r). For simplicity, we ignore the angular dependence 
in the following which only cause correction factors of the order of one [51]. Within 
the donor-acceptor model the exciton dipole moment is described by /i = \fl\ = fio\i){i\ 
where we furthermore have assumed that the dipole moments at each chromophore site 
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Figure 1. Donor and acceptor with two environmental field with different wave length 
are illustrated. The donor holds an electron in the excited state which forms with the 
hole left in the ground state an exciton having a finite electric dipole moment. 



are the same (again neglecting factors of the order one due to angular dependencies). 
The electric field is proportional to the amplitude of the propagating normal modes of 
the medium, E{r) = sgn{E{r)} ■ \E{r)\ = i?o(l/V^) ^kC*'''" finally leading to the 
interaction Hamiltonian in Eq. (|2]) with Ak(r) = {^qEq/ \/N)e^^^ . Similar ideas have 
recently been used to understand the phonon influence on double quantum dot charge 
qubits [371 EH] or tunneling defects [51] . 

For propagating modes in three spatial dimensions the spectral function of intersite 
fluctuations between chromophores i and j becomes 

Tf \ ^k(ri)A_k(rj) . . _ . sin(wto) 

k 

with to = — (9) 

V 

with the sound velocity v (assuming linear dispersion tUk = vk and k = |k|), the distance 
rij = \ri — Tjl between site i and j, coupling strength a and upper cut-off ojc using an 
exponential form for the cut-off function. For the small cut-off frequencies typical for 
biomolecular environments [uc — A) the cut-off function will modify quantitatively but 
not qualitatively the results. However, no detailed information about the specific cut- 
off functions for biomolecular environments is available in the literature. The (on-site) 
spectrum is Ohmic (linear in u) for the FMO complex [521 [23] . Linear dispersion for 
the normal modes is a strong assumption and it is not clear if the simple Debye picture 
holds in biological soft matter up to energies A. For on-site fluctuations at site i the 
spectral function simplifies 

M^) = E ^l^'^l^ - ^k) = 2au;e--/- . (10) 

Alternatively one might consider localized environmental fluctuations with 
localization length ^, which is taken to be independent of the mode energy u, resulting 
in the spectral density 

Jl]'{uj) = 2au;e-"/"=e-"'^/« (11) 
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Figure 2. Occupation difference of donor and acceptor versus time and speed of sound 
for Wc — A = 106 cm^^, T = 152 K, a — 0.08 and r^a — 3.8 A assuming propagating 
environmental modes. 



for the intersite spectrum whereas the on-site spectrum, Eq. (fTUj) . is unaltered. 

The intersite spectrum J\'j^{u) vanishes for sites far apart r^j ^ ^ and the 
fluctuations at the donor and the acceptor sites are uncorrelated. The intersite spectrum 
becomes identical to the on-site spectrum for close sites r^j ^ ^. In the later case, the 
environmental fluctuations of both sites are fully correlated and thus actually identical. 
As illustrated in Fig. [1] by the blue line, both excitons then "see" the same electric field, 
which, in turn, only modifies the total energy but not the energy difference between 
donor and acceptor. Accordingly, these fully correlated fluctuations cannot influence 
the dynamics of the donor-acceptor system. 

For the case of propagating modes, qualitatively the same holds. When Ato ^ 1 all 
modes with u; > A will not contribute to the intersite spectral function. At least at weak 
coupling (a <^ 1), mainly the modes resonant with the tunneling splitting are relevant 
and thus we expect the intersite spectrum to be irrelevant. The environment acts as two 
independent baths at each chromophore site, as usually assumed in the literature. When 
the shortest wave length Ac in the spectrum is larger than the distance between donor 
and acceptor, we have that ujcto = {rij/Xc) ^ 1. Then, the environmental fluctuations 
are fully correlated between the sites and thus do not affect the donor-acceptor system. 

For the donor-acceptor system we can transform the Hamiltonian as discussed 
in Eq. (|3]) and then define a single effective spectral function (as for a spin-boson 
problem [201 HH]) resulting in 



which is of super-Ohmic form. Super-Ohmic environmental fluctuations can neither 
cause overdamping (except at large temperatures) nor localization, in clear qualitative 
contrast to pure Ohmic fluctuations. This is a drastic qualitative effect which spatially 
correlated environmental fluctuations cause on coherent exciton transfer. 




(12) 



In the limit oOcto ^ 1, we get 




(13) 
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2.4- Dynamics of a single transfer step 

Having determined the reduced density matrix by QUAPI, we can evaluate the 
occupation difference P(t) = (cTz) between donor and acceptor. In Fig. [21 P(t) is plotted 
over time versus sound velocity v with which the modes are assumed to propagate. We 
have chosen all parameter to match rather closely the properties of chromophores in the 
FMO complex [52]. We have used A = 106 cm~^ as tunneling element which corresponds 
to the largest coupling in the FMO complex [52] between chromophore 1 and 2 but have 
neglected the energy difference between the two sites. The distance between site 1 and 
2 in the FMO complex of Chlorobium tepidum [7] is r 12 = 3.8 A which are the closest 
two chromophores. Site 2 and 7 are maximally apart, = 11.3 A. The bath cut-off 
frequency varies in the literature [521 153] between Uc = 32 cm~^ and 150 cm~^. To be 
specific, we choose Uc = 106 cm~^ and temperature T = 152 K= A/kB- 

We find that quantum coherent oscillations occur which decay within about 1 ps 
for the smallest value of the sound velocity of several hundred m/s. The environmental 
fluctuations are uncorrelated in this case. For larger sound velocities, meaning increasing 
correlations of the fluctuations, the decay slows down considerably as expected since the 
wave length of the modes causing decoherence becomes larger than the distance between 
the chromophores and thus cannot harm coherence any longer. We are not aware of 
experimental data regarding the precise values of sound velocities for the biological 
embedding materials of the FMO complexes. As a guide we might use the sound velocity 
of water, v ~ 1500 m/s, which falls into the range of our plot. When frozen to ice, as 
in the low temperature experiments at 77 K or 180 K, one finds that v ~ 3150 m/s and 
coherence lives considerably longer. 

In order to elucidate the dependence of this effect on the system-bath coupling. 
Fig. [3] a) shows the result for a = 0.2, resulting in a reorganization energy A ~ 2aoJc = 
42.5 cm~^. This compares to the case shown in Fig. [2l where we have set a = 0.08, 
resulting in a reorganization energy A ~ 2aujc = 17 cm"^. The stronger coupling results 
in faster decoherence and the crossover to fully correlated environmental fluctuations 
causing long-time coherence is pushed to larger values of the sound velocity. 

To address the temperature dependence, we show in Fig. [3] b) and c) the results 
for P{t) for the same parameters as in Fig. [21 except that temperature is set to 
T = 76 K in b) and T = 304 K in c). At lower temperatures coherence expectedly 
survives longer whereas at higher temperatures coherence lives shorter. Nevertheless, 
for all temperatures the profound effect due to the finite sound velocity is present. 
This turnover between uncorrelated (strong decoherence) and strongly correlated (weak 
decoherence) environmental fluctuations thus is only weakly dependent on temperature. 
Next, we discuss the dependence of the crossover on the localization length ^. Fig. [H 
shows P{t) for localized modes with localization lengths between 0-2 nm. Again as 
expected, for small localization lengths, the fluctuations at each site are uncorrelated 
and the occupation difference decays in less than a picosecond. Assuming a distance 
between donor and acceptor of = 3.8 A, the localization lengths plotted in 
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Figure 3. Occupation difference of donor and acceptor versus time and speed of 
sound assuming propagating environmental modes. Parameters are chosen to be 
=^ A = 106 cm-i and r^a = 3.8 A, and in a) T = 152 K, a = 0.2, in b) T = 76 K, 
a = 0.08, and in c) T = 304 K, a = 0.08. 




Figure 4. Occupation difference of donor and acceptor versus time and localization 
length ^ of localized environmental modes. Parameters are chosen as uJc = A. = 
106 cm-\ T = 152 K, a = 0.08 and rda = 3.8 A. 

Fig. m reach up to about four times the donor-acceptor distance. Then sizable 
correlations are expected and coherent oscillations for more than 2 picoseconds occur. 
In conclusion, assuming locahzed or propagating modes resuhs qualitatively in the 
same behavior. The used parameters are all taken for the case of the FMO complex 
and thus our results strongly suggest that spatial correlations of the environmental 
fluctuations due to finite propagation time of the modes strongly influence the decay 
of coherence in exciton transfer processes. In order to judge on the quantitative effect, 
a comprehensive experimental investigation of the environmental modes is needed, in 
particular, whether the modes are propagating or localized, and accordingly, whether 
the sound velocity and/or the localization length is relevant. So far, we have discussed 
the influence of spatially correlated environmental fluctuations on the coherence of a 
single donor-acceptor pair. These results can in principle be extended to a chain of 
more chromophoric sites without changing the physical picture qualitatively. 
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Figure 5. Probability Pdi.d2{t) of both exciton being at the respective donor 
sites versus time and sound velocity assuming propagating environmental modes. 
Parameters are chosen to be = A = 106 cm~^, T = 15.2 K, and a = 0.04. The 
distances are for a) r^a = 3.8 A and ?- = 38 A , and for b) r^a = 38 A and r = 3.8 A. 

3. Transfer in two chromophore pairs 

3.1. Model 

The FMO complex consists of three identical subunits, each of which consists of seven 
chromophoric sites and acts as a conductor for the excitons. Most likely, this structure 
has been optimized with respect to efficiency and seems to contain also some redundancy, 
which might be a measure of reliability in nature. In any case, the complex structure 
gives raise to the question whether a crosstalk between the subunits exists, and this even 
in a quantum coherent manner. In turn, the question whether this serves any purpose 
for functionality of increased efficiency is reasonable. 

In order to approach this question on a qualitative level, simple low-dimensional 
effective models are necessary. We model a single subunit by one donor-acceptor pair 
and discuss in the following two such donor-acceptor pairs located at a distance r. We 
assume that each pair initially contains a single exciton at the respective donor site. 
We explicitly suppress exciton transfer from one pair to the other and start from the 
Hamiltonian 

1 ^ ^ 
H^^^ = -AY,{\d^j){a,\+\a,){dj\} K) ^k(ri)gk 

j=l j=l i=aj/dj k 

+ + • (14) 

k 

We assume in the following that each chromophore couples separately to the 
environmental fluctuations but in contrast to the previous section, here there are two 
distinct distances involved. As before, each donor is separated by a distance r^a = 
l^ai —'"'dil = |ra2 "T^j | from its acceptor. We assume this distance to be identical for both 
donor-acceptor pairs. The separation between both pairs is r = Ir^i — FojI = Ir^i — rrf^l- 
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3.2. Results 

Since we have seen that propagating and locahzed (with a finite locahzation length) 
modes cause similar results for the decay of coherence, we restrict the following 
investigation to propagating modes. As before we set the tunneling element A = 
106 cm~^ and the fiuctuation cut-off frequency Uc = 106 cm~^. Fig. |5] shows the 
probability Pdi,d2{^) ^^at both excitons are located at the donor sites versus time for 
various values of the sound velocity for a weak coupling a = 0.04. Note that we start 
from Pdi^ii^ = 0) = 1. The reorganization energy is then A ~ 2aujc = 8.5 cm~^ and 
a rather low temperature T = 15.2 K is chosen. In Fig. [5] a), we fix the distances 
'T'da = 3.8 A and r = 38 A which refiects two rather distant donor-acceptor pairs. 
Both donor-acceptor pairs should thus be independent of each other since no direct 
coupling is assumed. We find a qualitatively similar behavior as for the case of a single 
pair as described in Section |2l Differences arise due to the weaker coupling and lower 
temperature. In contrast to the former case, we show the results for two close-by donor- 
acceptor pairs, r = 3.8 A but with large donor-acceptor distance rda = 38 A in Fig. [5] 
b). Here, we find that the effect of decoherence increases with increasing sound velocity, 
in clear contrast to the previous discussion in Section O Hence, increasing spatial 
correlations between the two donor-acceptor pairs destroy quantum coherence. 

For a quantitative investigation of this observation, we determine the associated 
decoherence rate F by fitting an exponentially damped cosine to the data for Pdi,d2i't) 
shown in Fig. |5] and plot it in Fig. [6] versus the sound velocity for several distance ratios. 
As fundamental distance scale, we use ro = 3.8 A, the distance of chromophore 1 and 2 
in the FMO complex of Chlorobium tepidum [7\. At the sound velocity v = 7600 m/s, 
a mode can travel this distance within the time of = 50 fs. 

When both donor-acceptor pairs are far apart, r = 10 Tq, and the distance between 
donor and acceptor is rda = (shown by the blue up-triangle in Fig. [6]), we recover 
the result of the previous Section |2l The decoherence rate decreases with increasing 
sound velocity. When the donor-acceptor distance is also large, rda = 10 rg, there is 
no dependence of F on f for the investigated range of sound velocities (yellow left 
triangles in Fig. [6]). The spatial correlations of the fiuctuations simply do not extend 
from the donor to the acceptor site and thus the fiuctuations are uncorrelated. When 
the two pairs are close to each other and the donor and acceptor sites are also close, 
the decoherence rate F decreases with increasing sound velocity (black circles in Fig. |6]). 
A totally different case is reached when donor and acceptor sites are well separated, 
rda = 10 To and rda = 100 tq (red squares and green diamonds in Fig. [6]), but the two 
pairs are close, r = tq. With increasing sound velocity, the decoherence rate increases 
and approximately doubles. The effect is slightly larger with larger distance between 
donor and acceptor. 

With increasing sound velocity, the wave lengths of the propagating modes increase 
and thus spatial correlations of these modes reach further. When the fluctuations at 
the sites of both donor-acceptor pairs are fully correlated, which is approximately the 
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Figure 6. Decoherence rate F associated to tlie occupation probability Pdi.d2 (t) versus 
sound velocity v. Parameters are chosen as A = 106 cm^^, T — 15.2 K, a = 0.04 and 
u!c — 106 cm""'^. Inset: T versus the ratio rda/i^o for v ~ 4560 m/s in the range 
ro < < 2ro. 

case for sound velocities of f ~ 7600 m/s for r = tq, then the fluctuations of the second 
donor (acceptor) equally influence the first one. Thus, the influence of the fluctuations 
effectively doubles which is reflected in doubling the decoherence rate. When donor 
and acceptors are closer, they as well become more correlated with increasing sound 
velocity, which, as discussed in the Section |21 results in decreasing decoherence rates. 
Thus, two effects are competing here. This also explains why for r^a = 10 r^ (red squares 
in Fig. [6]) the increase of the rate F is weaker than for r^a = 100 tq (green diamonds in 
Fig. E]). When both distances are equal, r^a = r = ro (black circles), the suppression of 
decoherence due to spatial correlations between donor and acceptor sites is the dominant 
effect, but it is weakened in comparison to the case of r^a = and r^a = 10 ro (blue up 
triangles). Hence, we find two regimes, in which either one of the two effects dominates. 
How sharp the crossover between the two regimes is, becomes visible when looking at 
the data for r^a = 2ro and r = tq. Here again, the decoherence rate F is increased by 
spatial correlations. Similar behavior (not shown) is also found for a smaller fluctuation 
cut-off frequency Uc = 53 cm~^ and otherwise identical parameters. The inset of Fig. [6] 
shows the change of the decoherence rate F when changing the donor-acceptor distance 
rda from ro to 2ro for a fixed sound velocity v = 4560 m/s. We find that in this regime, 
the change is linear in r^a- 

Quahtatively the same happens for very large fluctuation cut-off frequencies Uc = 
1060 cm^^ (as shown in Fig. [7^ for otherwise identical parameters). Note that in 
this regime, non-Markovian effects are suppressed. The increase (decrease) of F for 
correlated donor-acceptor pairs (for correlated fluctuations at donor and acceptor site) 
becomes sharper. The influence of temperature is studied in Fig. Ub), which shows the 
decoherence rate F for a higher temperature, T = 152 K, but otherwise for the same 
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Figure 7. Decoherence rate F associated to tlie occupation probability Pdi .d2 (t) versus 
sound velocity v. Parameters are A — 106 cm^^, a = 0.04, and in a) T — 15.2 K, 
ujc = 1060 cm-\ and in b) T = 152 K, uj^ = 106 cm-\ 



parameters as in Fig. |6l For small donor-acceptor distance, r^a = '"o? we find as before 
that the decoherence rate decreases with increasing sound velocity. However, the picture 
changes for large donor-acceptor distance, r^a = lOro- Then, the effect of decoherence 
is almost independent of the sound velocity, irrespective of the distance r between the 
donor-acceptor pairs. 

4. Discussion and conclusions 

To summarize, we have investigated the effect of spatially correlated environmental 
fluctuations on the quantum coherent transfer dynamics of excitons in donor-acceptor 
systems. Technically, spatially correlated environmental fluctuations can be included 
by adopting the numerical quasiadiabatic propagator path integral scheme. The 
fluctuations at each chromophore site generate a separate term in the Feynman- Vernon 
influence phase. The spatial correlations generate additional terms which then describe 
bath modes at different positions. Nevertheless, only the Feynman- Vernon influence 
phase is modifled and the general QUAPI approach still remains feasible. 

When the spatial correlations of the environmental fluctuations extend from the 
donor to the acceptor site, their decohering influence is strongly reduced since then, the 
energies at both sites are identically modifled. The energy difference between donor and 
acceptor is not changed but only the global reference energy fluctuates which does not 
influence the dynamics. The spatial correlations of propagating modes are characterized 
by their wave length which itself is determined by the sound velocity assuming linear 
dispersion. In contrast, the spatial correlations of localized modes are determined 
by their localization length. Although these two cases are in principle different, the 
qualitative effect on quantum coherence is the same and depends only on the ratio of 
distance between donor and acceptor and the correlation length given either by the 
localization length or by the wave length of resonant modes. 
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Two donor-acceptor systems in close proximity in addition show an increase in 
the decoherence rate in dependence on their spatial distance. When the donor and 
acceptor sites are far apart and their fluctuations are only weakly correlated, correlations 
of close-by donor-acceptor pairs become relevant when the distance between the two 
donor-acceptor systems becomes small enough that their environmental fluctuations are 
correlated. Then, each donor-acceptor system "sees" the fluctuations at the site of the 
other one and thus the decoherence rate is doubled. For intermediate spatial distances, 
both effects are competing with each other. The effect of an increased decoherence rate, 
however, is suppressed at higher temperatures and is probably less relevant at room 
temperature. 

Our results show that quantum coherence in the excitation transfer crucially 
depends on spatial correlations in the environmental fluctuations as soon as correlations 
lengths are of the order of the spatial distances of the chromophores. We have shown that 
for realistic material parameters (in particular, sound velocities), noticeable influence 
of the finite propagation times of environmental modes occurs. In fact, correlated 
fluctuations can reduce the decohering effect of the chromophore environment, when the 
correlation range extends over typical transfer distances within a chromophore chain. 
This effect even survives (even though diminished) up to room temperature and thus 
might be relevant for exciton transfer in biological systems. These findings are in line 
with the experimental findings [H |9]. 

At the same time, correlations between different chromophore chains increase the 
decoherence rate and thus a close packing would be disadvantageous. This increase is, 
however, suppressed at room temperature. Thus, increased thermal fluctuations actually 
indirectly support quantum coherence since they reduce spatial correlations which would 
increase decoherence. This effect could be tested experimentally by performing the 
reported electronic coherence photon echo experiments at even lower temperature at 
a full FMO complex with all three subunits. At low temperature, spatial correlations 
between the FMO subunits should reduce quantum coherence more effectively than at 
higher temperatures. 
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